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ABSTRACT 


;  An  asymptotic  method  is  developed  for  the  study  of  gravity-capillary 
waves  in  a  compressible  viscous  fluid  with  edge  constraints  in  an  inclined, 
straight  channel.  The  Navier-Stokes  equations  subject  to  free  surface  and 
rigid  bottom  conditions  are  reduced  to  a  sequence  of  elliptic  boundary 
problems  over  a  cross  section  of  the  channel.  Their  solutions  are  used  to 
determine  the  wave  speed  and  to  construct  the  Burgers  equation  for  the 
evolution  of  the  gravity-capillary  waves.  The  Burgers  equation  may  become 
ill-posed  when  the  Reynolds  number  exceeds  some  critical  value.  A  criterion 
for  the  stability  of  the  flow  is  then  defined  in  terms  of  the  critical 
Reynolds  number.  '  ' 
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SIGNIFICANCE  AND  EXPLANATION 


We  consider  a  compressible  viscous  flow  down  an  inclined  straight  channel 
of  arbitrary  cross  section.  The  fluid  is  filled  up  to  the  brim  of  the  channel 
so  that  the  edges  of  the  free  surface  of  the  fluid  remain  fixed.  Assume  that 
a  disturbance  of  small  amplitude  and  long  wave  length  is  present  on  the  free 
surface;  we  would  like  to  study  the  subsequent  development  of  the  disturbance. 
Within  the  framework  of  the  so-called  long  wave  approximation,  we  scale 
various  variables  by  appropriate  units  and  expand  the  solution  of  the  Navier- 
Stokes  equations  governing  the  motion  of  a  compressible  viscous  fluid  in  an 
asymptotic  series  in  terms  of  a  small  parameter.  The  governing  equations 
subject  to  prescribed  boundary  conditions  and  edge  constraints  are  reduced  to 
a  sequence  of  two-dimensional  linear  elliptic  boundary  problems  over  a  cross 
section  of  the  channel.  We  make  use  of  their  solutions  to  determine  the  wave 
speed  and  to  construct  the  Burgers  equation  for  the  evolution  of  the 
disturbance.  The  reduction  of  the  Navier-Stokes  equations  subject  to  various 
boundary  conditions  to  a  single  nonlinear  equation,  for  which  the  solution 
method  is  well  known,  is  our  main  contribution.  When  the  Reynolds  number  of 
our  problem  exceeds  some  critical  value,  the  Burgers  equation  may  not  be  well- 
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* 

M.  C.  Shen 


1.  Introduction 

In  recent  years  there  has  been  growing  interest  in  the  study  of  the  initial-boundary 
value  problem  of  a  compressible  viscous  fluid.  A  bibliography  up  to  1960  may  be  found  in 
the  review  article  by  Solonnikov  and  Kazhikov  (1981).  The  problem  of  a  compressible 
viscous  fluid  with  heat  conduction  bounded  by  a  rigid  boundary  was  resolved  by  Matsumura 
and  Nishida  (1981),  and  that  subject  to  free  surface  conditions  only  was  considered  by  Tani 
(1981).  However,  it  appears  that  the  general  case  with  both  free  surface  and  rigid 
boundary  conditions  still  remains  a  difficult  problem  and  it  may  become  more  untractable  if 
the  free  surface  and  the  rigid  boundary  intersect.  Needless  to  say,  at  present  no 
quantitative  analytical  results  are  available  for  these  problems.  On  the  other  hand,  the 
study  of  gravity-capillary  waves  on  an  inviscid  fluid  with  edge  constraints  in  a  channel 
has  received  some  attention  lately  (Benjamin  and  Scott,  1979;  Benjamin,  1981;  Shen  1982). 

It  should  be  of  great  interest,  therefore,  to  develop  an  asymptotic  method  for  the 
investigation  of  the  effects  of  coiqpressibility ,  viscosity  and  surface  tension  on  the 
propagation  of  surface  waves  subject  to  edge  constraints  in  a  channel.  The  problem  we 
shall  consider  deals  with  a  compressible  viscous  flow  down  an  inclined  straight  channel  of 
arbitrary  cross  section.  We  keep  the  edges  of  the  free  surface  fixed  by  filling  the  fluid 
up  to  the  brim  of  the  channel.  Suppose  a  surface  disturbance  of  small  amplitutde  and  long 
wave  length  is  present  in  the  channel.  We  would  like  to  study  the  subsequent  development 
of  the  disturbance. 
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Our  method  of  approach  is  based  upon  the  so-called  long  wave  approximation  (Shih  and 
Shen,  1975).  We  first  introduce  a  small  parameter  as  the  ratio  of  the  vertical  length 
scale  to  the  horizontal  length  scale.  Then  we  use  it  to  stretch  various  variables  in  the 
governinq  equations.  A  coordinate  system  is  chosen  to  move  with  the  disturbance  at  a  speed 
near  some  critical  speed  to  be  determined.  By  expanding  the  solution  of  the  governing 
equations  in  an  asymptotic  series/  we  reduce  the  Navier-Stokes  equations  with  boundary 
conditions  to  a  sequence  of  two-dimensional  linear  elliptic  boundary  value  problems  over  a 
cross  section  of  the  channel.  Their  solutions  are  then  used  to  determine  the  critical 
speed  and  the  Burgers  equation  as  an  approximate  equation  for  the  evolution  of  the 
disturbance.  The  Burgers  equation  may  become  ill-posed  if  the  Reynolds  number  defined  in 
our  problem  exceeds  some  critical  value.  We  make  use  of  the  critical  Reynolds  number  to 
define  a  stability  criterion  of  the  disturbance.  The  linear  stability  theory  of  an 
incompressible  viscous  flow  down  an  inclined  plane  was  first  studied  by  Yih  (1963),  and  the 
ill-posed  problem  of  the  linearized  Burgers  equation  for  the  same  flow  was  considered  by 
Carasso  and  Shen  (1977). 

We  formulate  our  problem  in  Section  1.  The  critical  speed  is  determined  in  Section 
2.  In  Section  3,  the  Burgers  equation  is  derived  and  the  stability  criterion,  defined. 

Some  special  cases  such  as  a  compressible  viscous  flow  down  an  inclined  tube  and  two- 
dimensional  flow  down  an  inclined  plane  or  between  two  planes  are  discussed  and  some 
remarks  regarding  a  heat -conductive  fluid  and  problems  for  further  study  are  given  in 


Section  4 


2.  Fornulation 


We  consider  the  notion  of  s  compressible/  viscous  fluid  in  an  inclined  straight 
channel  of  arbitrary  cross  section.  A  rectangular  coordinate  system  (x^x^x-j)  moving 
with  a  constant  speed  X*  in  the  x^-direction  is  chosen  so  that  the  two  edges  of  the 
channel  lie  in  the  x^/Xj-plane  and  H*(x2,Xj)  is  the  equation  of  the  channel  surface 
(Figure  1).  For  simplicity,  we  assume  that  the  fluid  is  barotropic,  that  is,  the  density 
p*  of  the  fluid  is  a  given  function  of  the  pressure  p*  only.  In  reference  to  the  moving 
frame,  the  Navier-Stokes  equations  are 

3p»/3t*  +  7*«p*q*  -  0  ,  (1) 

p*  Dq*/Dt*  =  -7*p*  +  p*q  +  u  7*2^*  +  (u/3 )7*(7*-q* )  ,  (2> 

P*  -  P*(p*)  (3) 


subject  to  the  boundary  conditions 
at  the  free  surface 


at  the  channel  wall 


at  the  edges 


£  *  -x3  +  nMXj.Xj.t*)  ”0  , 

°ijnj  “  T(Rt-1  +  Rf1)ni  ' 

D£/Dt*  -  0  f 


H*(x2,x3>  -  0  , 


q  -  <-X*,0,0)  , 

*2  -  tt*,  x3  »  0  , 

n*  -  o  . 


(4) 

(5> 


(6) 


(7) 


Here  7*  *  (3/3x^,  3/3x 
vector,  g  =  (g  sin  fl, 
inclination  angle  and 
stress  tewsor  given  by 


n  =  (3fi*/3x1,3nV3x  ,  - 


2<  3/3x3>,  D/Dt*  *  3/3t*  ♦  q»7*  q*  “  (q^q^qj)  is  the  velocity 
0,  -g  cos  9)  is  the  constant  gravitational  acceleration,  0  is  the 
0  <  0  <  x/2,  u  is  the  constant  viscosity  coefficient,  is  the 

=  (-p*  -  2y7*»q*/3)«ii  tuUq^x^  +  2qj/3xi>  ,  (8) 

1)  is  the  normal  vector  to  the  free  surface  £  =  0,  T  is  the 


constant  surface  tension  coefficient,  and 


4 


M 


n 


i 


»*-1  ♦  B*"1  -  tn*  d  ♦  n*2)  +  n*  <1  +  n*2)  - 

1  *  x2  *2*2  *  1 

2  2  -3/2 

2b*  n*  n*  ]  (i  ♦  n*  ♦  n*  )  , 

X1X2  X1  X2  X1  X2 

where  the  subscript  denotes  the  partial  differentiation  with  respect  to  . 

Furthermore,  we  also  assume  that  p*(p*)  >  0,  dp*/dp*  >  0. 

We  introduce  the  following  nondimensional  variables 

2  V 

t  =  e  t*/(h/g)  2 ,  (x,y,z)  *  fexf  ,x2 ,x3  )/h  , 

—  —i  V 

q  «  (u ,  v,  w )  «  (q  ^ ,  e  q2«  e  q3>/(gh)  2,  p  *  P*/A  , 


(9) 


p  =  p*/<Agh),  n  *  n*/b,  H(y,2)  »  H*(y*,2*)  , 

R  -  A(gh)1/2h/u,  S  =  T/Iu(gh>1/2  ],  X  »  X*/(gh>1/2  , 


t  -  i*/h,  e  -  h/L  «  1  , 

where  h  and  L  are  respectively  the  vertical  and  horizontal  length  scales,  A  is  a 

reference  value  for  the  density  to  be  chosen  later.  In  terms  of  them,  the  governing 

equations  and  boundary  conditions  (1)  to  (7)  become 

ep  +  (pu)  ♦  (pv)  +  (pw)  «  0  , 

t  x  y  2 

pe<eu  +  uu  ♦  vu  ♦  vu  )  *  -€p  +  p  sin  8 

t  x  y  z  x 

+  H~\e2u  +  u  ♦  u  )  *  R_1(e2/3)(V*q) 

XX  yy  22  x 

2  -12 
pc  (v  +  UV  +  UV  +  WV  )  a  -py  +  R  C  (C  v  +  V  +  V  ) 

t  x  y  y  —  —  — 


(10) 


Mi) 


xx  yy  zz 

+  R-1(C/3)(7-q> 

y 

2  -12 
pe  (w^  +uw  +  vw  ♦  ww  )  =  -p  -p  cos  0  +  R  e  (e  w  +  w  +w  ) 
t  x  y  z  z  xx  yy  zz 

+  R-1(e/3)(7»q) 

2 

p  =  p (p)  , 

subject  to  the  boundary  conditions:  At  the  free  surface  z  »  n(x,y,t,e). 


-4- 


(12) 


(13) 


(14) 


u 


k 


ti 


* 


?! 


t 


e  [Bp  +  (2  e/3)v»q  -  2eu  In  -  (u  ♦  e  v  )n 

x  x  y  x  y 

2  -1-1 
+  u  +  e  w  =  -s  (  r  +  r,  )en  , 

z  x  1  2  x 

[Rp  +  (2  c/3)7«q  -  2ev  1  n  -  e(e^v  +  u  )r\ 
y  y  x  y  x 

+  e ( v  +  w  )  =  -sir'1  ♦  r'1)!!  , 

z  y  1  2  y 

Rp  ♦  (2  e/3)7«q  -  2ew  +  elc^w  +  u  )n  +  e  (w  +  v  )n 

Z  X  z  X  y  z  y 

-1  -1 

-  -siP1  +  r2  )  , 

where  7«q  =  (u  ♦  v  +  w  ), 
x  y  z 

-1-12  2  22,2  ,  222  -3/2 

R,  +R-  *  en  ( 1  ♦  n  )  ♦  n  d  +  m  )  -  2m  nnd+en  +  n  ) 

12  xx  yyy  x  xy  x  y  xy 

en  +  un  +  vn  -w  =  o  ; 
t  x  y 


at  H (y , z )  »  0, 


at  y  =  ±1,  z  -  0, 


u  =  -X,v=w=0  ; 


n  =  o  . 


3.  Critical  speed 

Assume  that  the  solutions  q,  P,  p,  0  and  A  possess  an  asymptotic  expansion  of  the 


*  "  *0  *  £*1  +  6  *2  * 


We  call  Ag  the  critical  speed.  Without  loss  of  generality,  assume  A  ~  Aq  +  eA^. 
Substitution  of  (21)  in  (10)  to  (20)  will  yield  a  sequence  of  equations  and  boundary 
conditions  for  the  successive  approximations.  The  equations  for  the  zeroth  approximation 


are  the  following: 


-12 

R  7  Ug  «  -pQ  sin  0  in  r  , 
Uqz  »  0  on  L0  , 
u_  =  -A.  on  L .  , 


We  obtain  from  (36)  and  (37)  that 


p1z/p1  "  ■O'(P0)  coa  0  • 

By  integration  and  making  use  of  (25),  we  have 

P,  “  P(p„>f ^x.t)  , 

where  fj(x,t),  a  function  of  x,  t  only,  is  to  be  determined.  Next  from  (25),  (39) 
(42)  and  (43)  it  follows  that 

Sn«  -  R  cos  Bn,  “  -Rf,(x,t) 

lyy  1  1 

n1  =  0  at  y  =  ±Jt 

It  is  easy  to  obtain  that 

n  1  =  v,)  (y)f  ^x.t)  , 

where 


v 


1 


(y >  =  (cos  9  )  1  [1 
k  »  (R  cos  9/S)^2 


-  (cosh  kl )  cosh  ky) 

2  1/_ 

“  (pgh  cos  9/T)  2  . 


Let 


u,  ■  -  X1  • 


Then  from  (35),  (38)  and  (41),  4^  is  the  solution  of  the  following  problem: 


V  4,  *•  'D*(p0>p0  sin  9  in  D  , 


4.  *  ♦„  vj  -  4.  v,  on  L  , 

1*  Oy  1  T0zz  1  0 


4  ^  •  0  on  L 1  . 


To  determine  XQ,  we  integrate  (34)  over  D,  use  the  divergence  theorem,  (30),  (37), 


(40),  (41),  (43)  and  (44)  to  obtain 

'l0[*o  -  Wix*  ■  -V[po*i  +  p’(po>po<Rpo  ■  Vlf1x“ 

and 


X0  ■  p,(p0)P0d*  +  /L  vi<Syl_1  t/ c/  <po*  1  +  p'(p0>P0*0>dA  +  ^ L  Vo*1 

0  0 


which  yields  the  expression  for  the  critical  speed.  We  write 


X.  -  RX*  , 
0  0 


where  X*  Is  independent  of  R. 


4.  Burgers  equations 


Our  next  step  is  to  derive  the  Burgers  equation  for  f.j(x,t),  and  we  proceed  to  the 
equations  for  the  second  approximation: 


P1t  +  <P0U2  +  P1U1  +  P2U0)x  +  <P0V2  +  P1V1>y  +  (P0W2  +  P1W1 


)  =  0 
z 


P0(U0U1x  *  V1U0y  +  W1U0z)  =  _P1x  *P!Sin!tR  7  U2 
p2y  =  r'VVi  +  ( 1/3  )R_  1  (V*  q  )y  , 


p2z  =  -p2  cos  9  +  R  V2w^  +  (1/3)R  ^(7»q2> 


2 '  z.  ' 


at  z  *»  0, 


at  H(y,z)  »  0, 


at  y  a  tl,  z  =  0, 


p 2  =  p  '  (P0>P2  +  P“<P0>P/2 


n„  -  u«  *1 «  +  u.  n4/2  +  u,.  t)  +  u-  =o 

Oy  2y  1y  1y  Ozz  2  Ozz  1  Izz  1  2z 

-  uA  n«  +  v,  +  w  *  o  , 

0y  lx  1z  1y 

Rp2  *  RP0zn2  +  Rpizni  +  (2'3>7^1  '  "  -SR2yy  ' 

nit  +  Uinix  +  V2X  +  Viniy  =  W2  +  W1zni  f 


u2  =  v2  =  v2  =  0  ’ 
n2  “  0  ' 


We  first  go  back  to  (34)  and  make  use  of  (37),  (43)  and  (46)  to  express  (34)  as 

<poVy  +  <P0W1)z  =  "Vi  +plVx  (63) 

-  -R[poS  +  (*o  -x*o)p,,po>po]f,x 

Then  we  construct  a  function  Q(y,z)  as  a  solution  of  the  following  equations: 

(p qVq,)  -  -p^,  +  <♦„  _xo)p,(po>po  in  °  '  (64> 

°1Z  ‘  <*„  “  XS>V1  °"  L0  '  (65> 

Pn  =  0  on  .  (66) 

The  Neumann  problem  posed  by  (64)  to  (66)  is  solvable  because  of  (50).  From  (63),  (64),  we 


<P0V1  •  P0RC1yf1x)y  +  !P0W1  "  Wz'lx*  “  °  • 

It  follows  from  (67)  that  we  may  define  a  function  R4.(x,y)  such  that 


(67) 


I 


R4,  f,  “  P„V,  -  pRO  f  , 
1z  lx  0  1  0  1y  lx 

-  R4,  f,  »  P„W4  -  p  RQ  f  . 

1y  lx  01  0  1z  lx 


Hence , 


V1  "  R(p0  *12  +  °1y)f1x  ' 


1 


R(-p_1*  +  Q  )f 

1y  1z  lx 


(68) 

(69) 


We  now  cross-differentiate  (54)  and  (55)  and  subtract  one  equation  from  another  to  obtain 

-1  2 

R  7  (v,  -w,  )  =  -p_  cos  9 

I2  1y  2y 

By  substituting  (68)  and  (69)  for  and  w1  in  the  above  equation  and  making  use  of 

(46),  (54),  (68)  and  (69),  it  is  obtained  that 


V2  (V*  (pQ1*1  ))  =  -cos  0p'  (pQ  )  [V2  (pQ1*1i )  +  72Q1y  +  ( 1/3  )*  1 
+  +  (po’*1y,Z  +  72q1]  in  °  * 


(70) 


The  boundary  conditions  for  (70)  can  be  derived  from  (40),  (41),  and  (58).  It  follows  from 
(30),  (44),  (68)  and  (69)  that 


-1 


iy 


=  0,  <p0  *u)z  =  ♦  vt  -  2[(*0  -  X*)^]y  on  LQ  , 


P„  +  Q, 

0  1z  1y 


°'  Po’*1y 


°1*  “  °  L1 


(70) 


Suppose  Q1  and  4 ^  have  been  found,  and  in  the  following  we  shall  successively 
construct  manageable  expressions  for  p2>  hj  and  Uj.  We  multiply  both  (54)  and  (55)  by 
PQ  and  make  use  of  (56)  to  obtain 


<p01p2>y  “  P01b  1  +  n/3)  (7*q1  )y]  , 


(72) 


(pg’pj'jj  =  Pg1  t-p"(pQ)p2  cos  0/2  +  R_172Wi  +  ( 1/3 )R_  1  (7» )z J  •  (73) 

In  D,  a  fixed  point  (x,y0,Zg)  and  a  smooth  curve  V  from  (Xg,y0,Zg)  to  <Xg,y,z) 
are  chosen.  Then  from  (72)  and  (73),  we  obtain 

p2  =  Pq  / 2  PolR  Vv,  ♦  <1/3)(7*q1)y]dy 

+  P-1 t-p"(p. )p?  COS  0/2  +  R  '72w  +  (1/3)R  1(7»q)z)dz  +  p  f  (x,t) 

u  U  1  I  u  01 


"  P20f1  +  P21f1x  +  Pnf 


0  2 


(74) 


where  f2  is  an  unknown  function  of  x  and  t,  and  by  (43),  (46),  (68)  and  (69), 
p20  “  ■(p0/2>  ^ L  P"<P0)P0  cos  0  ' 

P21  *  p0  U  pi1{v2tp0l41«  +  V  +  (V3>f*ly  +  <p.’,*l.  +  01yV 


+  +  *  {^(Pfl^ly  +  ®1z>  + 

+  '-*•**  +  <V«*,}d2  • 

Next  we  turn  to  the  equation  (59)  for  hj.  By  (25),  (36),  (37),  (43),  (44),  (46) 
(62),  (68),  and  (69),  we  may  express  (59)  as 


2yy 


(R/S)  cos  Mj*  (R/S)[-p2Q  +  p'(0)  cos  8 1  f  1 

<R/S){p21  4-  (2/3)  (f ,  ♦  (p;1*^  *  01y>y]  -  *  ClB)B)f1x 


n2  «  0  at  y  “  tt 


2  2  2 

Let  G(y,z)  be  the  Green's  function  for  (d  /dy  -  k  )n2  with  n2  “  0  at  y  ”  tt , 


easily  found  that 


Then 


G(y,z)  ■  sinh  k(C-t)sin  k(y+£)/(k  sinh  2kt),  -t  <  y  <  C  , 
sinh  klc+t  )sinh  k(y-t)/(k  sinh  2ki),  C  <  y  <  l  . 


where 


ho  “  v,.f,  +  v..f,  v,f,  , 
2  20  1  21  tx  12' 


v2o  =  /-t  G(1-')WS)<T>20  tp'(0)  cos  0lz»o«iy  > 

V21  =  !-l  G<y*«)<R/S>fP21  +  <2/3>  ,  ♦  tpQ1*!*  +  Ply>yJ 


-  (d/SK-p;1*^  ♦  0, z>z>z-0dV  • 

Finally  we  take  up  u2.  Making  use  of  the  expressions  derived  before  for  u0,  q^ 
p2,  ni  and  n  ,  we  obtain  from  (53),  (56),  (57)  and  (61)  that 


it  is 


(75) 


i  P)  * 


-11- 


R  7  u2  -  {R  P0I<V*5>*1  +  <p0  *1*  +  W  +  ^0  %  +  P1*>80*J 

+  P0  "  P21  “in  8  P'(p0^f1x  "  [p’(p0>p20  +  P"  (P0)P0/2]«inc  9  f2 
-  Pg  sin  9  f2  in  D  , 

U2*  “  R(80yV21  ”  ♦o*«V21)f1x  +  W80yW20  +  ♦lyVi  '  80**V20 

"  *0*x*V?/2  *  ♦l2ZV1,f2  *  R(*0yVJ  -  80**Vf2  °n  S  ' 

u_  •  0  on  L. 

2  1 

As  observed  from  the  above  equations,  we  may  express  u2  in  terms  of  f1x,  f  and  f2 
follows: 

u2  “  [R3923(y,*>  +  "♦21<y'*,Jf1x  +  R*3<y**)*i  ♦  R*.,f2  •  ( 

Here  *2#  *3>  *4  are  the  solutions  of  the  following  problems: 

V2*23  “  P0I(80  “  X0,81  +  (p0l81*  +  °1y^0y  +  ^O^ly  +  Qls^Ox1'  in  D  ' 

8 23*  “  °  °n  L0  ' 

*23  -  0  on  L,  , 

72*2i  -  P0  -  P21p,(P0)  «in  9  in  O  , 


6  m  a  v  •  -  A  v  onL  * 

V21s  v0y  21  *0*s  21  1  ' 


*21  -  0  \  » 


V*3  -  -[p,(p0)p20  +  p"(p„)P0/2]sin  9  in  D  , 

8 3*  “  80yv20  +  81yV1  ”  80**V20  ~  80***V1/2  “  81**V1  on  L0 


*3*0  on  L1 


Now  we  are  in  a  position  to  derive  the  Burgers  equation  for  f As  before  for  the 


derivation  of  the  critical  speed,  we  integrate  (52)  over  D  and  apply  the  divergence 
theorem  and  (61)  to  obtain 


IjJ  [(e0V2,y  +  “W.1®  *  k0  W2dY 

(71 

"  ~U  P1t  *  (P0U2  +  Pl“t  *  P2U0,*dA  '  \  P1W1dy  ' 

where  we  note  that  pQ  •  1  on  Lq.  By  making  use  of  (30),  (37),  (40),  (43),  (44),  (51), 
(56),  (60),  (74),  (75)  and  (76),  and  rearranging  the  terms,  (77)  becomes 

B0f1t  *  “lf1x  +  m2ff1x  ”  m3f1xx  '  (7e 

Here  the  coefficient  of  f2  vanishes  because  of  (50),  and 

mo  “  U  v,dy  + //D  p,(p0)podA  '  179 

0 

mi  *  imQ  '  <8C 

”2  ■  R/L0t2(VX0,W20  *  V,,po’*1*  +  Ply’  •  U1(“p01*1y  + 

,,  (81 

♦  p,<p0HVX0,V1ldy  *  ffD(2p0*3  +  *"(P0,P0*1 

+  2p,(p0)(*0-x5)P2o  ♦  <>“<Po,pO(VX0,,dA  ' 

”3  "  [p0,r3*23  +  *21*  *  Rp’(p0)p21(VX0)ldA 

-  /l0W*0**0>W21^  • 


(78)  with  (79)  to  (82)  is  our  main  result.  It  is  well  known  that  the  Burgers  equation  is 

not  well  posed  if  m3  becomes  negative.  We  set  m3  “  0  and  solve  for  R  to  obtain 

Rc  “  critical  Reynolds  number 

-1  *  1/ 

-  (I -JDJ  f>o*23dAl  [V  lp0*21  *  P’<p0,p21(VXO,)dA  +  /L0(pO-X0,V21dyl}  2* 

Assume  Rc  is  positive  and  finite  and  m3  <  0  when  R  >  Rc.  Then  the  Burgers  equation  is 
ill-posed  for  R  >  R  and  we  define  R  <  R  as  a  criterion  for  the  stability  of  the  flow. 


5.  Discussion 


In  the  following  we  consider  sooe  special  cases,  which  can  be  easily  dealt  with  by  our 
previous  results. 

(1)  Fixed  boundary  case 

This  case  corresponds  to  a  compressible  viscous  flow  down  an  Inclined  tube  of  uniform 
cross  section.  He  only  have  the  boundary  condition  u  “  -X,  v  *  w  ”  0  on  h(y,z)  »  0.  In 
all  the  elliptical  boundary  value  problems,  we  extend  the  boundary  condition  on  L.,  to  the 
whole  boundary.  Evidently  in  (79)  to  (82)  all  the  line  Integrals  on  Lg  should  be 
dropped. 

(2)  Two-dimensional  case 

In  this  case,  the  solution  of  the  governing  equations  (1)  to  (6)  is  independent  of  the 
coordinate  y,  and  we  also  assume  v  =  0.  Furthermore,  the  effect  of  surface  tension  will 
not  manifest  itself  up  to  the  equations  for  the  second  approximation.  First  we  consider  a 
compressible  viscous  fluid  flow  down  an  inclined  plane.  p0  is  still  determined  by  (29) 


and  pfl,  by  (28).  (31)  to  (33)  become 


+  “  “  -pQ  sin  0  in  -1  <  z  <  0  , 


♦o 


♦o  -° 


at  *  *  o  , 
at  z  »  -1  , 

where  0g  is  a  function  of  z  only.  For  the  first  approximation,  pj  is  given  by 
(43).  (44)  should  be  replaced  by 


h  j  m  f  ^  sec  0 

Noting  that  =  0  in  (38),  (47)  to  (49)  now  become 

♦  ^  *  -p'(pg)pg  sin  0  in  -1  <  z  <  0 


(83) 


♦j  “  *0  8ec  9 

*1  ■  ° 


at  z  "  0  , 
at  z  “  -1 


He 


integrate  (34)  without  (p^v^)^  respect  to  z  from  -1  to  z  to  obtain 

W1  ’  ■p()1^-1<P0u1  +  plVxdl 


(84) 


FF(z)f 


lx 
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where 


*•<*>  ■  -Oq1  /!i[po*i  +  p,(po,po{*o  '  xoR’1]<5* 

At  z  »  0,  (40)  holds  and  It  follows  from  (83)  and  (84)  that 

X.  -  RX9  , 

0  0 


X*  -  (/°t  p’(p0)P0dz  +  sec  8)"1t/®1(p0*1  +  p  •  (p„  )P040  )dz  +  ♦{)(0)sec  9]  . 

Since  the  expression  for  w1  is  known,  we  need  not  consider  the  equations  (63)  to 
(71).  From  (74)  we  obtain 

p2  =  p0/!l  p01f*p"(p0,p1  008  0/2  +  R"1w1zz  +  (1/3,R"1(u1x  +  w1z,zld* 

+  “oV*'*’  ( 

"  P20f1  *  P21f1x  +  P0f2 


p20  ’  "<P0/2)  ^-1  P“(P0,P0  COS  9  dZ  ' 
p2 1  -P0  Jf,  P;’tM/3)f  +  (1/3)9']dz  . 


With  similar  changes  in  (59),  we  have 


n2  *W20f1  +  V21f1x  +  (8,C  ®,f2  ' 


V20  "  tp20<0)  "  P*<0,,88C  8  * 

V21  “  tp21(<,>  +  <2/3)9,(0)  -  (4/3 )F’  (0)  ]sec  8  . 

Finally,  as  in  the  three  dimensional  case,  we  express  U2 

u2  "  *  *21(*,,f1X  +  "♦3<l,f1  +  "*1*2  ' 


*23  -p0t(VX0’*1  *  in  -’  <  *  <  1 

*23  “  0  at  1  "  0  ' 


0  at  z  “  -1  1 


*21  *  P0  "  p21p,<p0>Sin  6  in  -1  <  *  <  1  ' 

*21  "  "*0^21  at  *  ’  °  ' 

*21  “  °  at  *  “  _1  ' 

♦  J  -  ~tP'(P0)P2o  *  <*  ”*P0  JP  q/2]  Bin  8  in  -1  <  z  <  0  , 

*3  ”  "*^20  “  *0  s«c2®/2  ‘  MC  8  at  *  *  0  * 

■  0  at  z  ■  -1 

To  derive  the  Burgers  equation,  we  integrate  (52)  without  (pQv2  +  pivi*y  wittl  respect 
to  z  from  -1  to  0  and  obtain 

m0f1t  +  m1f1x  +  m2f1f1x  “  n3f1xx  ' 


mQ  -  sec  8  ♦  J®,  p'(p0»P0dz 
«  sec  8  +  (p  0  ( — 1 )  -  1]sec  8 
“  P0(-1)sec  8  , 

"l  *  "  xiB'o  ' 

■j  -  2f(80(0)  -  XJ]u20  -  rf* (0)sec  8 

+  /!1I2p0B*3  *  2R(P--Xo,p,(p0,p20  *  *’<VP0  +  p"(p0>p0,d*  ' 

m3  -  -R[V0)  -  X*]v21  -/!1tP0(l«3823  ♦  *21>  ♦  BP*(p0)(80-X5)p21]dz 
As  before,  we  set  *3  -  0  *nd  solve  R  to  obtain 

Rc  "  tt"/!lp0*23a*1’1^-1to0*21  +  p,{p0,p21<*0_X0)]da 

+  <V°>  -X5,V211)1/2  • 


Assume  Rc  is  positive  and  finite.  If  for  R  >  Rc,  n>3  is  negative,  then  we  define  R  < 
Rc  as  the  criterion  for  the  stability  of  the  flow. 

Suppose  the  plane  z  "  0  is  a  rigid  boundary.  Then  in  all  the  two-point  boundary 
value  problems,  we  simply  impose  at  z  ■  -1  the  same  boundary  condition  as  at  z  ■  0. 
The  coefficients  in  the  Burgers  equation  become 
m_  •  tp.l-l)  -  1]  sec  8  , 


»2  ’  /!,f20*3  *  2B(VX0,P'<P0)P20  +  ^,<p0)p0  *  p“<p0)p0ld*  ' 

"3  “  "f-1  tp0(1,3423  *  "♦21)  +  rp'(p0)(VX0)p21ldl  • 

In  conclusion,  we  make  a  few  remarks  regarding  some  possible  extension  of  the  method 
developed  here  and  problems  for  further  study.  If  we  assume  the  fluid  is  heat-conductive, 
then  the  relation  p  »  ptp)  should  be  replaced  by  the  energy  equation  involving 
temperature  and  appropriate  boundary  conditions  for  teiqperature  should  be  prescribed  at  the 
free  surface  and  rigid  boundaries.  The  same  method  can  be  carried  through  if  the 
coefficient  of  heat  conductivity  is  sufficiently  large.  The  critical  case  R  «  Rc 
certainly  is  an  interesting  problem  and  it  should  be  of  importance  to  derive  an  asymptotic 
equation  near  and  at  R  ■  Rc  to  replace  the  Burgers  equation.  Furthermore,  if  R  >  Rc, 
the  Burgers  equation  may  be  ill-posed  and  a  study  of  the  ill-posed  problem  is  also  of 
significance. 
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